Atomic states in optical traps near a planar surface 
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In this work we discuss the atomic states in a vertical optical lattice in proximity of a surface. 
We study the modifications to the ordinary Wannier-Stark states in presence of a surface and we 
characterize the energy shifts produced by the Casimir-Polder interaction between atom and mirror. 
In this context, we introduce an effective model describing the finite size of the atom in order to 
regularize the energy corrections. In addition, the modifications to the energy levels due to a 
hypothetical non-Newtonian gravitational potential as well as their experimental observability are 
investigated. 
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I. INTRODUCTION 



Atomic interferometry has the potential to become 
a powerful method to investigate atom-surface interac- 
tions, the main reason being the high precision vifhich can 
be reached in frequency measurements. In this context a 
new experiment named FORCA-G (FORce de CAsimir 
et Gravitation a courte distance) has been recently pro- 
posed [l|. The purpose of this experiment is manifold: 
on one hand it aims at providing a new observation of 
the Casimir-Polder interaction between an atom and a 
surface, resulting from the coupUng of the fluctuating 
quantum electromagnetic field with the atom on the 
other hand it also intends to impose new constraints on 
the existence of hypothetical deviations from the New- 
tonian law of gravitation. These goals will be achieved 
thanks to the innovative design of FORCA-G, in which 
interferometric techniques are combined with a trapping 
potential. This is generated by a vertical standing optical 
wave produced by the reflection of a laser on a mirror. 
The vertical configuration leads to an external potential 
on the atom given by the sum of the optical one and a 
linear gravitational term due to the earth: this deviation 
from a purely periodical potential produces a localization 
of the atomic wavepacket, corresponding to the transi- 
tion from Bloch to Wannier-Stark states The main 
advantages of FORCA-G are thus the refined control of 
the atomic position as well as the high precision of in- 
terferometric measurements, as demonstrated in the first 
experimental results [1]. 

Having in mind a theory-experiment comparison 
within a given accuracy, the theoretical treatment of the 
problem as well as the experimental investigation must 
be independently assessed with the same precision. In 
the case of FORCA-G this demands a detailed theoreti- 
cal study of the atomic wavefunctions and energy levels 
in proximity of a surface. As an intermediate step, a pre- 
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cise characterization of the Casimir-Polder atom-surface 
interaction is also needed. 

These issues are the main subject of investigation of 
this paper. As a matter of fact, the influence of the 
Casimir-Polder interaction on the atomic energy levels 
has so far been explored [ll, Q using the simple idea of 
calculating the electrodynamical potential at the center 
of each well of the trap. We will discuss the validity of this 
model focusing in particular on the scheme of FORCA- 
G. In this work, we present a hamiltonian approach to 
this problem. This treatment flrst allows us to discuss, 
independently on the Casimir-Polder atom-surface inter- 
action, the atomic trapped states. Since the presence of 
the surface breaks the translational symmetry typical of 
Bloch and Wannier-Stark problems, we focus in partic- 
ular on the difference (both in energy levels and wave- 
functions) between our trapped states and the standard 
Wannier-Stark solutions. Then, in order to discuss the 
Casimir-Polder corrections to the energy levels, we gener- 
alize the perturbative treatment usually exploited to de- 
duce atom-surface electrodynamical interactions, by in- 
cluding the external optical and gravitational potentials, 
and treating as a consequence the atomic coordinate as 
a dynamic variable. The theoretical work presented here 
will be useful for all experiments that aim at measuring 
short range interactions between atoms trapped in opti- 
cal lattices and a macroscopic surface [1, 5, 6] as they will 
require precise modelling of the atomic states and energy 
levels close to the surface. 

This paper is organized as follows. In section |TT] we 
describe our physical system. Then, in section HIIl we dis- 
cuss the shape of atomic wavefunctions in the trap. Sec- 
tion |IV] is dedicated to the study of the Casimir-Polder 
interaction and its influence on the atomic energy lev- 
els. In this section we introduce an effective description 
of the finite size of the atom and discuss its validity in 
connection with the experiment. In section |V] we look 
at the energy shifts introduced by a hypothetical non- 
Newtonian potential and we investigate the constraints 
that FORCA-G could impose on the strength of this de- 
viation. Finally, in section IVll we discuss our results. 
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II. THE PHYSICAL SYSTEM 

In this section we are going to describe the main fea- 
tures of our physical system and the hamiltonian formal- 
ism used to investigate the interaction between atom and 
electromagnetic field. Let us consider a two-level atom 
trapped in an optical standing wave produced by the re- 
flection of a laser having wavelength Xi — ^ on a. surface 
located at 2: = 0. In the configuration we are consider- 
ing the optical trap has a vertical orientation, so that 
we have to take into account the earth's gravitation field 
acting on the atom. The complete hamiltonian can be 
written under the form 

H = Ho + Hint = Hf + i/at + Hws + Hint 

p + OO p 

Hf = Y^ J dk^ J cfkhujal(k,k^)ap(k,k^) 



H^t = hhJo\e){e\ 

p2 U , 



(1) 



^^ws = ^ + "^52 + "^(1 ~ cos(2/c/z)) 

Hint = ■ f(r)- 

The complete Hamiltonian is written as a sum of a term 
Hq describing the free evolution of the atomic and field 
degrees of freedom. In particular, iJf is the Hamiltonian 
of the quantum electromagnetic field, described by a set 
of modes (p, k, fc^): here p is the polarization index, tak- 
ing the values p = 1, 2 corresponding to TE and TM po- 
larization respectively, while k and kz are the transverse 
and longitudinal components of the wavevector. We as- 
sociate to each single mode a frequency w = c-\/P + AJf , 
as well as annihilation and a creation operators ap(k, k^) 
and aj^(k, fcz). An eigenstate of the field Hamiltonian is 
thus specified by giving a set of photon occupation num- 
bers |{np(k, kz)}) for each mode of the field. The vacuum 
state of the field, with zero photons in each mode, will be 
noted with |Op(k, fcz)). In our formalism the expression 
of the electric field is the following 



• /.+00 
£{v) = -^J dkz j d^k 
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(2) 



X e 



ikrj 



fp(k, kz,z)ap{\^, kz) - h.c.) 



where we have introduced the transverse coordinate rj^ = 
(x,?/) and the mode functions fp(k, fc^jz) characterizing 
the boundary conditions imposed on the field. Under 
the assumption of a perfectly conducting mirror in z = 
these functions take a very simple expression ^] 

fi(k, kz,z) —\<iX isn\{kzz) 

. (3) 

f2(k, kz, z) — k sm{kzz) — z — cos(fczz) 

UJ UJ 

where k = k/fc and z — (0,0,1). H^t is the inter- 
nal Hamiltonian of our two level atom having ground 
state \g) and excited state |e) separated by a transition 



frequency uiq- While i?at is associated to the internal 
atomic degrees of freedom, the term TJws accounts for 
the external atomic dynamics. As a consequence, it con- 
tains the kinetic energy {p being the canonical momen- 
tum associated to z), as well as both the gravitational 
potential (treated here in first approximation as a linear 
term), where m is the atomic mass and g is the accel- 
eration of the Earth's gravity, and the classical descrip- 
tion of the stationary optical trap, having depth U. We 
treat here only the z-dependent terms of the Hamiltonian 
since the degrees of freedom x and y, even in presence 
of a transverse trapping mechanism, are decoupled from 
the longitudinal dynamics. For simplicity, we shall take 
as a unit of energy the photon recoil energy given 

by Er — -2^- As far as the atomic position is con- 
cerned, it will be expressed in units of the periodicity 
of the trap ^■. For all numerical examples in this pa- 
per we will use the experimental configuration chosen for 
FORCA-G: Er = 5.37 x lO^^o j ^ g.ll x 10^ Hz) 
and ^ = 266 nm. 

The interaction between the atom and the quantum 
electromagnetic field is written here in the well-known 
multipolar coupling in dipole approximation Q, where 
fi — qp (q being the electron's charge and p the internal 
atomic coordinate) is the quantum operator associated to 
the atomic electric dipole moment and the electric field 
is calculated in the atomic position r. It is important to 
observe that, since /x clearly operates only on the atomic 
internal states, this interaction term is the only one cou- 
pling atomic (both internal and external) and field de- 
grees of freedom. As a consequence, the ground state of 
the free Hamiltonian Hq is simply given by the tensor 
product of the vacuum field state |Op(k, kz)), the atomic 
state \g) and the ground state of Hy/s- In the picture of 
atomic dressing [9|, the ground state of Hq is the bare 
ground state, and the inclusion of Hint will produce a 
new ground state of the complete system, referred to as 
dressed ground state, mixing all the degrees of freedom. 
We are going to tackle the calculation of the ground state 
of -ffvvs in the next section, whereas the atom-field inter- 
action will be treated in section HVl 



III. MODIFIED WANNIER-STARK STATES 

A. Ordinary Wannier-Stark states 

In solid state physics, it is well-known that the solution 
of the time-independent Schrodinger equation describing 
a quantum particle in a periodic potential leads to the 
so-called Bloch states [ifl, Due to the periodicity 

of the system, these states are completely delocalized in 
space coordinate and the spectrum energy is composed 
of bands of permitted energies, each band being labeled 
with an index b = 1,2,.... The addition of a linear po- 
tential (whose role is in our case played by gravity) to the 
trap produces localization of the states: these states are 
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usually labeled as Wannier-Stark states (see e.g. [I2|,|l3j) . 
We will now describe their main features. For each Bloch 
band b, a discrete quantum number n is introduced, tak- 
ing the values n = 0, ±1, ±2, The state \n, 6)*^^' is, 

in coordinate representation, approximately centered in 
the rt-th well of the optical trap, and the energy of this 
state is in first approximation given by 



1^^0,1 (z) 



p(WS) 



p(WS) 



mgXi 



(4) 



with i;^"^' the average of the 6-th Bloch band 0, 
As a result of the quasi-periodicity of the system (i.e. 
of the linearity of the gravitational potential modifying 
the periodic trap) two states jn, fe)*™^' and |p, 5)'™^' be- 
longing to the same band h are shifted, in coordinate 
representation, by n — p wells. At the same time their 
energies differ, in accordance with eq. (|4]), by n — p times 
5g — . Then, the problem of Wannier-Stark states is 
solved once we know, for each band b, the average Bloch- 
band energy E^^' and the eigenfunction centered in a 
given well. The Wannier-Stark states can be calculated 
using, for example, the numerical approach of [isj . 

In figures [T] and [2] we give the Wannier-Stark states 
|0, 1)'^^' for two different values of the potential depth 
U = 3,10 (in units of Er). The figures show that, as 
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FIG. 1: Coordinate representation of the state |0, 1) belonging 
to the first Bloch band and centered in the zeroth well for 
U = 3. 

expected, a deeper well produces a more localized state 
of the particle. 



B. Wannier-Stark states in proximity of a surface 

In the context of our problem, the presence of a surface 
at z = plays two roles. On one hand, it induces a 
modification of the Wannier-Stark states by imposing a 
boundary condition on the eigenvalue problem. On the 
other hand, the quantum electrodynamical interaction 
between the atom and this surface must be taken into 
account, as we will describe in section ITVl 

The surface aX z = breaks the quasi-periodicity of the 
system. The potential modifying the optical trap is no 
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FIG. 2: Coordinate representation of the state jO, 1) belonging 
to the first Bloch band and centered in the zeroth well for 
U = 10. 



longer linear, since it must be considered as the gravita- 
tional linear potential for z > and an infinite potential 
barrier for z < 0, describing the impossibility of the par- 
ticle to penetrate into the mirror. We will refer to the 
eigenstates of this new physical system as the modified 
Wannier-Stark states. From now on we are going to deal 
only with these new states: the state of the 6-th Bloch 
band centered in the n-th well will be noted with \n,b) 
(and correspondingly ■0„^b(z)). 

We have solved the problem of modified Wannier-Stark 
states numerically, using a finite-difference method. The 
first step of our approach consists in considering a unidi- 
mensional box < z < Zf and imposing that the wave- 
function vanishes at the borders. As for z = 0, this cor- 
responds to a real physical boundary condition, whereas 
the condition ip{zf) = is purely numerical. Naturally, 
the acceptability of the solutions will depend on their 
rate of decay toward for z — >■ z/. The next step is the 
discretization of the interval [0, z/] using a set of + 2 
mesh points Zi with zq = 0, zi, . . . , zat+i — Zf (giving 
Sz = for equally spaced mesh points). 

Using this approach, the problem is reduced to an 
eigenvalue problem of a tridiagonal symmetric matrix. 
The solution of such a problem can be efficiently worked 
out using the numerical approach first introduced in [isj 
as well as a standard QL algorithm [3]. In order to 
check the robustness of our numerical results, we have 
also checked their coherence with a finite-element method 

Choosing a large enough numerical box, taking for ex- 
ample z^ = 30 (we recall here that z is measured in units 
of trap periods ^), we have verified that the modified 
Wannier-Stark states centered in a well far from the sur- 
face (approximately starting from n = 10) have the same 
shape as the functions shown in section lTlI Al this refiects 
the fact that far from the surface the quasi-periodicity of 
the system is reestablished. Moreover, in this region we 
find that the energy difference between two successive 
states equals the expected quantity Sg defined before: 
starting from n = 10, the differences equal Sg with a rela- 
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tive precision better than 10~^. This can be seen in table 
|T1 where we show the resuhs obtained for the first ten en- 
ergy levels with [7 = 3: in this table we give the energy 



levels En, as well as the differences 5 En = En+i — En 
In this configuration we have 6g = 0.070068. 
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FIG. 3: (Color online) Density probability of modified Wannier-Stark states tp„^i{z) for n = 1,2,3,4 and [7 = 3. The last two 
functions (black, solid line) are compared to the corresponding standard Wannier-Stark state (red, dashed line). 



Table U shows only values of the energies belonging to 
the first Bloch band {b = 1). We have checked that, in- 
creasing the value of N above 10^, the last digit reported 
in table H] remains constant, corresponding to a relative 
precision of approximately 10"''. Actually, as a result of 
our numerical method we also found the eigenvalues and 
corresponding states associated to higher bands. Never- 
theless, we shall discuss only the first-band states since 
the ones belonging to higher bands are much less relevant 
for experimental purposes: as a matter of fact, higher 
bands are not efficiently trapped in the experiment (the 
average energy of the second band is around AEr for a 
trap depth of 3Er). 

As far as the states in proximity of the plate are con- 
cerned, they are strongly modified by the boundary con- 
dition, and the same property holds for their energies. 
We show, in figure |3] the first four eigenfunctions in pres- 
ence of the surface. For the sake of comparison, the third 



and fourth wavefunctions are superposed to the standard 
Wannier-Stark solutions centered in the corresponding 
well. It is important to stress that the ordinary Wannier- 
Stark functions of wells n = 3, 4 are plotted only to show 
that the shape of the modified ones tends towards the 
standard solution: however, the fact that the ordinary 
functions for these wells are different from zero for z < 
makes strictly speaking no sense for our physical system. 

In order to discuss the influence of the depth of the 
wells, we will conclude this section giving the results 
obtained for U = 10. In this case, since the ordinary 
Wannier-Stark states are much more localized in each 
well, we expect the influence of the surface to be evident 
on a smaller range of distances. This can be seen directly 
from table ini where the energy differences converge more 
rapidly to Sg. 

Moreover, from figure 0] we see that the state 4'2,i{z) 
shows already a remarkable accordance with the corre- 
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n 


En 


SE„{xlO ^) 


(5iJ„(xlO^Hz) 


1 


1.4028 


12.302 


9.9788 


2 


1.5258 


9.8043 


7.9525 


3 


1.6239 


8.4432 


6.8485 


4 


1.7083 


7.6206 


6.1812 


5 


1.7845 


7.2026 


5.8422 


6 


1.8566 


7.0518 


5.7199 


7 


1.9271 


7.0146 


5.6897 


8 


1.9972 


7.0079 


5.6843 


9 


2.0673 


7.0070 


5.6835 


10 


2.1374 


7.0068 


5.6834 


11 


2.2074 


7.0068 


5.6834 


12 


2.2775 


7.0068 


5.6834 


13 


2.3476 


7.0068 


5.6834 



TABLE I: First ten values of the modified Wannier-Stark 
spectrum for f/ = 3. These values have been obtained on 
an interval [0,30]. The first two columns are in units of Er, 
the third one is in Hz. 



n 


E„ 


SE„{xlO-'^) 


(5£'„(x102Hz) 


1 


2.9496 


7.5127 


6.0938 


2 


3.0247 


7.0276 


5.7003 


3 


3.0950 


7.0072 


5.6837 


4 


3.1651 


7.0068 


5.6834 


5 


3.2352 


7.0068 


5.6834 


6 


3.3052 


7.0068 


5.6834 


7 


3.3753 


7.0068 


5.6834 


8 


3.4454 


7.0068 


5.6834 


9 


3.5154 


7.0068 


5.6834 


10 


3.5855 


7.0068 


5.6834 



TABLE II: First ten values of the modified Wannier-Stark 
spectrum for U = 10. Same parameters as for table H] 



spending uninodified Wannier-Stark state in the interval 
[0.6,6], where the probability of finding the atom is ap- 
proximately 0.9997. 



IV. CASIMIR-POLDER INTERACTION 
A. Standard Casimir-Polder calculations 

The presence of the surface does not only play the role 
of imposing a boundary condition on the Wannier-Stark 
wavefunctions. In fact, since it modifies the structure 
of the modes of the quantum electromagnetic field, it is 
source of an attractive force between the atom and the 
plate. This is a particular case of a general phenomenon 
usually called Casimir effect for two macroscopic bodies 
and Casimir-Polder force when it involves one or more 
atoms near a surface (for a general review see e.g. [l9{). 
This phenomenon was first pointed out by Casimir in 
1948 for two parallel perfectly conducting plates [13] and 
in the same year by Casimir and Polder for atom-surface 
and atom-atom systems 21 1 . 

The Casimir-Polder force between an atom and a mir- 
ror has been measured quite recently using several dif- 
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FIG. 4: (Color online) Density probabihty of modified 
Wannier-Stark states '4}2.i{z) (black, dashed line) for f/ = 10 
compared to the corresponding standard Wannier-Stark state 
(red, dotted). 



ferent techniques: deflection of atomic beams [22], reflec- 
tion of cold atoms 23-25]. In the past few years, Bose- 
Einstein condensates proved to be efficient probes of this 
effect, both by means of reflection techniques [1^ [l^] 
and observing center-of-mass oscillations of the conden- 
sate [28'-'3T|. The FORCA-G experiment aims at achiev- 
ing a percent precision in the measurement of the force 
thanks to the combination of cold atoms and interfcro- 
metric techniques. 

From a theoretical point of view, the force is usually 
obtained from an interaction energy which results from a 
time- independent perturbative calculation on the matter- 
field hamiltonian interaction term 0, |^ . In this kind of 
approach, the position of the atom is usually treated as 
a fixed parameter and not as a quantum operator. As a 
consequence, in order to deduce the Casimir-Polder inter- 
action energy between an atom and a perfectly conduct- 
ing plate, we must neglect the term Hws in the Hamil- 
tonian of the system ([T]) and use eq. ^ for the electric 
field. Choosing the bare ground state |0p(k, fc^)) ](;) as 
the unperturbed configuration, the first-order perturba- 
tive correction on interaction term 7?int is zero, since the 
atomic electric dipole moment operator is an odd opera- 
tor and the annihilation and creation operators appearing 
in the electric field do not connect states with the same 
number of photons. Moving to second-order, we obtain 
the z-dependent potential energy 



E 



+ 00 



^,^ |AW(k,/c.,r)|^ ^ (5) 



H{oj + Wo) 



In this expression we have defined 

A«(k,fc„r) = (Op(k,fc,)l(glHi„t|lp(k,fc,))le) 



TT V 47re 







and we sum over all the possible intermediate states 
[lp(k, fcz))]e) having one photon in the mode (p, k, fc^) 
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and the atom in its excited internal state |e). Finally, 
the superscripts (2) and (1) refer to the order with re- 
spect to the electric charge contained in fi. 

This result holds for a perfectly conducting surface and 
at zero temperature. However, the generalization to more 
realistic configurations including the finite conductivity 
of the plate as well as a temperature T > is not straight- 
forward in a perturbative approach. This can be worked 
out using for example the scattering method [H, [s^ or 
the Green- function formalism (see [SJ, HH] and references 
therein) . The resulting potential can be put under the 
form 1361 



(2) , , 2kBT ^ 2 «(^Cn) 



n=0 



dk 

'2Kn 



rTE{k,i^n) - ( H T2—]'^™{k,i^n) 



(7) 



£2 



where £„ = is the n-th Matsubara frequency and 

the prime on the Matsubara sum indicates that the n = 
term is to be taken with half weight. Moreover we have 

defined Kn — \J ^ + k^ and the rp{k,uj) are the well- 
known Fresnel coefficients for a planar surface. Finally 
a{ijj) is the ground-state atomic polarizability, which for 
a multilevel atom takes the form [33 



should be added to the Wannier-Stark Hamiltonian i?ws 
in (m to obtain a new time-independent problem. So 
one could obtain a new set of energies and wavefunctions 
taking also into account the quantum electrodynamical 
part of the problem. In figure[5]we plot this new complete 
potential (sum of ^ and the ordinary Wannier-Stark 
potential) for a Rubidium atom in front of a perfectly 
conducting surface at zero temperature. This plot clearly 
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FIG. 5: (Color online) Sum of Wannier-Stark (for U = and 
Casimir-Polder potentials (in black, solid line) compared to 
Wannier-Stark potential alone (in red, dashed line). 



a(aj) 



2 \ ^ ^-raOMwO 



(8) 



where £",10 = En — Eq is the difference between the en- 
ergies of the n-th atomic level (starting from the first 
excited state) and of ground state, whereas /i„o is the ma- 
trix element of the electric dipole operator between the 
same couple of states. Clearly, the conductive properties 
of the surface material are included in the Fresnel coef- 
ficients through the electric permittivity and magnetic 
susceptibility e{uj) and ^(w) respectively. We conclude 
this section giving the expression of the Casimir-Polder 
potential for an atom in front of a real surface at zero 
temperature 



+00 



-t-00 



-2Kz 



fTE(fc,i£) 
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dk- 



2K 



rTM(fc,i£) 



(9) 



where K = ^ ^ + k'^ and the sum over the Matsubara 
frequencies is replaced by an integral. 

B. Perturbation of modified Wannier-Stark 

In the last section we have given the Casimir-Polder 
potential for an atom having polarizability q;(w) in front 
of an arbitrary planar surface and at temperature T . It 
could be natural to think that this z-dependent potential 



shows that the Casimir-Polder interaction modifies the 
optical trap on a limited range. In particular, in our case 
the first well does no longer exist, the second and the 
third are slightly modified, and starting from the fourth 
the trap is practically unperturbed. 

Nevertheless, the simple addition of the Casimir-Polder 
z-dependent potential to the external hamiltonian term 
-ffws is strictly speaking incorrect. As a matter of fact, 
the potential ([5]) (as well as ([7]) and dU) has been de- 
rived using several hypotheses. First, it arises from a 
perturbative treatment of the interaction term Tfjut on 
the Hamiltonian iJat + Hi. Moreover, in this calculation 
the atomic position z is treated as a fixed parameter. 
This is clearly incoherent with the fact that for our com- 
plete Hamiltonian ([1]) z is a dynamic variable. 

These arguments suggest that we should reconsider the 
calculation of the Casimir-Polder potential using our per- 
turbative approach including the Wannier-Stark Hamil- 
tonian i/ws ■ The perturbative term is still i?int , but now 
the atomic coordinate z has to be treated as a quantum 
operator as well. So, we are able to introduce a new un- 
perturbed state for each well n of the first Bloch band 
having the form 



|VS) = |0,(k,fc,))|5)Kl). 



(10) 



As we found for ordinary Casimir-Polder calculation, the 
leading-order correction to the energies is the second, and 
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the corrections takes the new form 



(2) 



E 



dk. 



p 

+00 +00 

p=i 5=1 



(VSl^int|lp(k,fc.))|e)b,6)| 



(11) 



E. 



(0) 



(0) 
ri.l 



h{uj + UJo) 



where now the intermediate state contains the modified 
Wannier-Stark state \p,b). We notice that the difference 
between two Wannier-Stark energies appearing in the de- 
nominator is p—n times approximately O.OTEr if the state 
\p,b) belongs to the first band {b = 1). We now point 
out that the recoil energy for a Rubidium atom having 
TO = 1.44 • 10~^^ kg trapped in a periodic potential hav- 
ing A/ = 532 nm is of the order of 10~^^ eV. On the other 
hand, the atomic transition energy fkuo is of the order of 
the eV. At the same time, the numerator in (fTTj) involves 
an integral over the z coordinate containing the product 
of the wavefunctions associated to the two states. This 
product becomes negligibly small for |p — 7i| ^7, whilst 

the energy difference E^^J — E^f^^ is still orders of magni- 
tude smaller than Hloq. As a consequence, the Wannier- 
Stark energy difference in the denominator can be always 
safely neglected with respect to hujo for the intermediate 
states having 6=1. As far as the higher bands are con- 
cerned, it is possible to see that the same superposition 
integral decays to zero due to the delocalization of the 
modified Wannier-Stark states. Furthermore, in the case 
of higher bands, the energy difference E^J^j - e';°\ is still 
lOOEr (and then still negligible with respect to huJo) for 
b = 10. This reasoning enables us to use the closure 
relation on the \n,b) states and obtain 



(n,l|(-^ dk.Jd^k- 



4^Hk,fc,,r)|' 

h{uj + OJq) 



\n,l). 
(12) 



The expression in parentheses coincides with the second- 
order perturbative calculation on the atom-field ground 
state described in section lTV Al It is thus evident that the 
correction we are looking for equals the average on the 

Wannier-Stark state \n, 1) of the known Casimir-Polder 

(2) 

potential V^^{z). This can be then expressed as follows 



r+00 

AElll = / dz 
Jo 



V^^iz). 



(13) 



This expression has been obtained in the context of a per- 
turbative treatment for a perfectly conducting surface at 
zero temperature. Nevertheless, the reasoning which led 
us from the general expression (jlip to the simple average 

value ([T^ does not depend on the details of the calcula- 

(2) 

tion of the V^p (z) itself. As a consequence, it is reason- 
able to assume that the average value (IT51) can be also 



used with the more general expressions of the interaction 
energy or ©. 

The behavior of the integrand function around z = 
must be treated with care: indeed, the Casimir-Polder 
potential diverges for z 0. In particular, it is well- 
known that for distances much smaller than the typical 
atomic transition wavelength (van der Waals regime) the 
interaction potential is temperature-independent and its 
expression reads 



(2) 



CP;vdW 



(^) 



Attz^ 



+°° a(»0 6(»g) - 1 

^ 47reo e{tO + 1 



(14) 



e{uj) being the electric permittivity of the surface mate- 
rial. As for the atomic wavefunction, we have numeri- 
cally verified that, for any allowed value of n, it tends to 
zero linearly for z — > 0. As a consequence, the integrand 
function behaves like z~^ around the origin, implying a 
divergent energy correction (I13p for any n. We will de- 
velop in the next section an effective description of the 
atom to regularize this quantity. 



C. Regularization of the correction 

The potential V^p (z) represents a particular case of 
singular potentials since it diverges around the origin 
faster than z~^. The treatment of such potentials has 
been discussed since the pioneering work of Case [111 (for 
more details see e.g. [30|)- It can be shown from first 
principles [iOl that these potentials describe an unphys- 
ical situation in proximity of the origin. The solution 
of the time-independent Schrodinger equation requires 
in these cases a more detailed knowledge of the short- 
distance physics of the problem. 

In the case of atom-surface interaction, the z~^ behav- 
ior of the potential is an artefact of treating the atom as 
a point-like source. This statement is supported by the 
calculation of the Casimir potential between a sphere of 
radius R and a wall pll - l43| . The potential energy as- 
sociated with this geometrical configuration shows a z~^ 
long-distance behavior (equivalent to the long-distance 
Casimir-Polder atom-surface interaction), an intermedi- 
ate z~^ regime and a transition toward a z~^ behavior 
when approaching z = 0. This property holds for any 
value of the radius R of the sphere: nevertheless, the 
characteristic distance at which the transition occurs is, 
as physically predictable, of the order of the radius. 

Inspired by fiS], we take into account the finite size 
of the atom by replacing it with a probability density 
distribution p(r'). In accordance with [44|, we make the 
further assumption that the function p(r') is different 
from zero within a finite volume. Moreover, in our nu- 
merical applications, we take this volume to be a sphere 
of radius R, discussing also the dependence of the re- 
sults on R. We assume that the atom has coordinate 
(0,0, z). We stress here that the atomic coordinate z is 
taken at the point of the sphere nearest to the surface: as 
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a consequence the effective sphere representing the atom 
is centered in (0, 0, z + R). As far as the probabihty den- 
sity distribution is concerned, we wiU consider the cases 
of a constant function pi(r') — Ni and of a spherically 
symmetric parabolic distribution 



a 



P2{r') R^- - - (z' ~ z - Rf 



(15) 



For both probability distributions the variable r' is ex- 
pressed in the same frame of reference as for the atomic 
coordinate. The factors A^i and A^2 are to be deduced 
from the normalization condition 



d^r' p{r') 



1 



(16) 



being fl the spherical atomic volume. Our hypothesis 
leads to a new regularized expression of the atom-surface 
potential, given by the average with respect to p(r') of 
the standard Casimir-Polder potential 



^CP;reg(^)= / dh' piv' )Vg^ {z' ) 



(17) 



where the z-dependence of the new potential is implic- 
itly contained in the probability density distribution p(r') 
and the integration volume ft. Substituting p7)) into 
then provides the regularized energies of our system. 

Let us now analyze the behavior of the regularized po- 
tential (jl7l) in proximity of the surface. Assuming that 
it has a form 



(2) 



CP;rcg 



(Z) = — 

z" 



the exponent a has the form 



(2) 



CP;ro! 



dz 



(2) 



CP;rcs 



(18) 



(19) 



In figure [S] we plot the exponent a as a function of z 
for the standard Casimir-Polder potential ([9]) and the 
regularized one ([T7| . Both are calculated in this case 
for a Rubidium atom in front of a perfectly conducting 
surface and at zero temperature: the data for the dy- 
namical atomic polarizability of Rubidium were kindly 
provided by Derevianko et al. 45]. Furthermore, the 
regularized expression is calculated for a uniform proba- 
bility density distribution and three different radii R — 
100 pm, 1 nm, 10 nm. In the four cases it is evident that 
the transition from z~^ to z~^ behavior starts around the 
first atomic transition wavelength (~ 780 nm). Moreover, 
while for the standard Casimir-Polder calculation the ex- 
ponent tends to 3, in all the other cases the finite size of 
the atom leads, as anticipated, to a z~^ asymptotic de- 
pendence. The figure shows clearly that the length scale 
of this second power-law transition is roughly of the order 
of the atomic size. We will make use of this regularized 
potential in the next section to work out the perturbative 
calculations on the modified Wannier-Stark states. 



0.1 



1 



10 



100 



1000 



z[nm] 



FIG. 6: (Color online) Exponent a defined in H19[) for the 
standard Casimir-Polder potential (black, solid line), and 
the regularized potential (|17() for radii 0.1 nm (red, dashed 
line), 1 nm (blue, dotted line) and 10 nm (orange, dotted- 
dashed line). 



D. Energy corrections 

We are now ready to evaluate the average value of 
the potential (|17p on any modified Wannier-Stark state. 
Our approach leaves as free parameters the atomic ef- 
fective radius R and the probability density distribution 
p(r'). As far as the radius R is concerned, we first re- 
mark that several non-equivalent definitions of the effec- 
tive atomic radius exist in literature. For example Slater 
gives in (46| an empirical value of Rubidium radius equal 
to 235 pm with an associated accuracy of 5 pm. On the 
contrary, the work [Jzj estimates the atomic radius for 
Rubidium to be 265 pm. As a consequence, in order to 
study the dependence of the results on the value of the ra- 
dius, we will consider the two extreme cases R = 200 pm 
and R — 300 pm. As for the probability distribution, 
we will use the functions pi(r') and P2(r') discussed be- 
fore. Furthermore we are going to consider the case of a 
perfect conductor for the surface in order to get an in- 
sight on the qualitative features of the energy correction. 
In table IIIII we show the energy corrections to the first 
twelve modified Wannier-Stark states obtained choosing 
two radii and two probability density distributions. In 
order to be coherent with the description of the atom 
as a sphere, the same regularization treatment used for 
the Casimir-Polder interaction is applied to the other z- 
dependent hamiltonian terms contained in -ffws- As ex- 
pected (i? <C Xi) this does not modify our results by more 
than 10~^ in relative value. 



It is easy to see from table IIIII that, as far as the en- 
ergy levels are concerned, a change in the effective radius 
produces a more remarkable effect than a change of dis- 
tribution from the uniform case pi(r') to the parabolic 
P2{r'). In particular, switching from 200 pm to 300 pm 
gives a relative error which is of the order of 10% on 
the first wells and then drops down, whereas the relative 
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n 


200 pm - pi 


200 pm - p2 


300 pm - pi 


300 pm - p2 


1 


2.39[1] 


2.37[1] 


2.19[1] 


2.16[1] 


2 


1.76[1] 


1.74[1] 


1.60[1] 


1.58[1] 


3 


1.20[1] 


1.18[1] 


1.09[1] 


1.08[1] 


4 


6.80 


6.7147 


6.18 


6.09 


5 


2.89 


2.8557 


2.63 


2.59 


6 


8.71[-1] 


8.60[-l] 


7.91[-1] 


7.80[-l] 


7 


1.93[-1] 


1.91[-1] 


1.76[-1] 


1.74[-1] 


8 


3.57[-2] 


3.53[-2] 


3.27[-2] 


3.23[-2] 


9 


6.76[-3] 


6.7Q[-3] 


6.33[-3] 


6.27[-3] 


10 


1.84 [-3] 


1.83[-3] 


1.78 [-3] 


1.78[-3] 


11 


8.36[-4] 


8. 35 [-4] 


8.29[-4] 


8. 29 [-4] 


12 


5.10[-4] 


5.1Q[-4] 


5.09[-4] 


5.09[-4] 



TABLE III: Absolute value (in Hz) of the Casimir-Polder en- 
ergy corrections (they are all changed in sign) to the first 
twelve modified Wannier-Stark states for [7 = 3. The nota- 
tion a [6] corresponds to a x 10*". These values are calculated 
for a perfectly conducting surface. 



correction from pi(r') to P2(r') is at most around 1%. 

It is now instructive to compare one of the set of energy 
corrections shown in table Hill with the simple evaluation 
of the strength of the potential energy (0) at the centre 
of each well [l|, which could be used as a first esti- 
mation of the energy correction. This idea works better 
considering deeper traps or farther from the surface: for 
example, we have verified that the value of VJjp (z) cal- 
culated aX z = 1 and the first energy correction differ 
by a factor of approximately 4.4 for [/ = 3 (see figure [7]), 
while this factors drops already to 1.12 for [/ = 20 and to 
1.05 for U = 80. We also remark that a larger value of U 
or a larger atom-surface distance reduces the dependence 
of the results on the choice of both the probability den- 
sity distribution and the effective radius. This reasoning 
proves that in the case of the experiment FORCA-G the 
delocalization of the atom indeed plays a role. 

AE;(Hz) 

10 Z ~ - _ 

I — 
O.I - ~ 

O.OI - _ _ 

0.001 ~ - Z - _ 

2 4 6 8 10 12 ^ 

FIG. 7: (Color online) Absolute value of the Casimir-Polder 
energy correction for a uniform distribution and a radius of 
200 pm (black upper ticks) compared to the Casimir-Polder 
potential by evaluating ((7]) at the well centre (red lower ticks). 
The depth of the trapping potential is (7 = 3. 

We want to stress here that the validity of our 
spherical-atom model used for the regularization of 



V^'^iz) still remains to be tested by experimental mea- 
surements. Some more details, as well as the relationship 
with the search for non-Newtonian gravity, will be given 
in section IVTl 



V. DEVIATIONS FROM NEWTONIAN 
GRAVITATION 

Many theories of unification of general relativity and 
quantum mechanics predict a modification of the laws of 
gravity at short distances. These modifications can be 
described by the addition of a new potential to the stan- 
dard Newtonian one. This correction is often modelized 
by a Yukawa-type law so that the complete gravitational 
potential between two point-like particles is written un- 
der the form 

, , GMm ( , , 

[7(3 (z) = — - — (^l + aye ^yJ (20) 

where G is the gravitational constant, m and M the 
masses of the two particles. In this expression ay and 
Ay are two parameters introduced to characterize respec- 
tively the relative strength of the corrective potential 
and its typical range. The experiments aimed at test- 
ing the existence of such a deviation set constraints on 
the allowed values of the parameters ay and Ay. The 
present status of the excluded regions at short ranges 
(z < 100 /im) on the (ay. Ay) plane is depicted in figure 

m 

In the experimental configuration of FORCA-G we 
have verified that the only relevant Yukawa-type con- 
tribution is the one associated to the atom-mirror grav- 
itational interaction. At the same time, the Newtonian 
part of the atom-surface interaction is completely negli- 
gible with respect to the Earth- atom term already taken 
into account in the Wannier-Stark Hamiltonian ([Ij and 
with respect to the expected experimental uncertainties. 
As a consequence, the correction we are looking for is 
obtained by integrating the Yukawa part of eq. (1^ over 
the volume occupied by the surface. Describing the mir- 
ror as a cylinder (the atom being on the direction of its 
axis) and recalling the we are looking for deviations hav- 
ing length scale Ay in the /im range we obtain, after a 
straightforward calculation, 

Hy{z) = 2TTaYGpsm\\e'^ (21) 

ps being the density of the surface. We are now going to 
find the new unperturbed energy levels of the system (in 
absence of Casimir-Polder interaction) in presence of the 
new hamiltonian term (j21l) . This can be done using the 
method described in section lllli after having chosen the 
value of the parameters ay. Ay and ps- The new eigen- 
values of the unperturbed Hamiltonian i?ws + Hy will 
be noted with for each well n. As far as the surface 
density is concerned, since we still do not have any infor- 
mation about the surface to be used in the experiment. 
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we choose throughout this section just as an example the 
density of sihcon ps — 2.33 x lO^kgm"^, close to the 
values corresponding to Si02 or BK7 typically used in 
experiments. 

As anticipated in the introduction, one of the the 
scopes of the experiment FORCA-G is to look for 
Yukawa-type deviations both near the surface (at dis- 
tances of the order of /ini) and in the region where the 
Casimir-Polder interaction can be theoretically modelled 
at a degree of precision comparable to the experimental 
noise. In the former regime, the idea of the experiment is 
to compare the results obtained using two different iso- 
topes of Rubidum (in particular, *^Rb and *^Rb) in or- 
der to make the energy differences between wells almost 
independent on the Casimir-Polder interaction [Tl]. As a 
consequence, when discussing the Yukawa correction near 
the surface, we first need to calculate (both for Wannier- 
Stark and Yukawa potentials) the differences in energy 

levels En and between ^^Rb and ^''Rb, calculated 

using the formalism described in the previous sections 
with the different isotope masses in the Hamiltonians ([1]) 
and (I21|l . These differences will be noted with. 

VEn = - Efj) - - £;(^)«^) (22) 

Finally the experiment will be able to detect a Yukawa- 
type deviation if the difference "DEn is within the exper- 
imental sensitivity. In the case of FORCA-G, the ex- 
pected sensitivity is 10~^Hz [l|. 

We first give in table |IV] the results obtained for ay = 
3 X 10^'^ and Ay = 1 /im: the value of ay approximately 
corresponds to the limit of the experimentally accessed 
region for Ay = 1 A*™- 



n 


VE„ (Hz) 




n 


VE„ (Hz) 


1 


1.654[-1] 




13 


2.0 [-3] 


2 


1.425[-1] 




14 


1.5 [-3] 


3 


1.221[-1] 




15 


1.2 [-3] 


4 


9.43[-2] 




16 


9 [-4] 


5 


5.72[-2] 




17 


7 [-4] 


6 


2.71 [-2] 




18 


5 [-4] 


7 


1.30[-2] 




19 


4 [-4] 


8 


8.0[-3] 




20 


3 [-4] 


9 


6.0[-3] 




21 


2 [-4] 


10 


4.4[-3] 




22 


2 [-4] 


11 


3.4[-3] 




23 


l[-4] 


12 


2.6[-3] 




24 


l[-4] 



TABLE IV: First 24 values (the energy differences are ex- 
pressed in Hz) of the modified Wannier-Stark spectrum for 
(7 = 3 and in presence of the Yukawa-type potential (|2ip 
with qy = 3 X 10^" and Ay ~ Ifim. 



From these results it is clear that the Yukawa-type devia- 
tions corresponding to the couple (ay. Ay) chosen are, in 
principle, experimentally detectable up to the well n = 24 
in a differential ^^Rb-^^Rb measurement. 



We now turn to the second experimental configura- 
tion in which the Casimir-Polder potential is expected 
to be predicted at the 10~^Hz level [l|. We stress that, 
apart from the precision in the calculation of the non- 
regularized Casimir-Polder potential, we must pay at- 
tention to the uncertainty introduced by our effective de- 
scription of the finite size of the atom. Assuming that 
the Casimir-Polder potential can be theoretically deter- 
mined, independently of its regularization, with at best a 
1% accuracy, the absolute precision in its determination 
can be considered comparable to the experimental error 
already around the well n — 40, where z = 10 /im and 
the potential equals approximately 0.06 Hz, i.e. in this 
second experimental configuration the atoms will be at 
10 fim or more from the surface. At this distance, our 
hypothesis of a spherical atom plays already no role. In- 
deed, we have checked that, using both probability distri- 
butions and both radii, the potentials so obtained differ 
less than 10^^ Hz already aX z — b /im. This is coher- 
ent with the fact that the finite size of the atom plays a 
negligible role at distances much larger than the atomic 
size itself. As a consequence, in this second experimen- 
tal regime the precision on the standard calculation and 
the experimental uncertainties impose stronger limita- 
tions than our effective model. 

We have calculated the Yukawa corrections on the well 
71 = 40 for different values of Ay: for each of them, 
we have found our limiting value of ay by looking for 
a correction of the order of 10~^Hz. We have, more- 
over, repeated the same calculation for n = 70 (where 
Vcp — 0.01 Hz) as well as in the near regime discussed 
above evaluating the energy difference between wells 
n = 4 and n = 6. These three curves are represented in 
figure [5] on top of the present experimental constraints. 



VI. DISCUSSION 

In this work, we have introduced an effective model de- 
scribing the atom as a spherical probability distribution. 
This was needed in order to regularize the expression 
of the Casimir-Polder energy correction to the modified 
Wannier-Stark states. It is important to discuss in more 
detail the validity of this model in connection with the 
experimental results. 

Let us start by recalling that in the context of the 
search for non-Newtonian gravitation, our model does 
not impose severe limitations. As a matter of fact, in 
the near regime the Casimir-Polder contribution is al- 
most cancelled by the use of two isotopes, whereas at 
far distances we have shown (see sec. |V]) that the error 
introduced by our description is negligible with respect 
to the accuracy in the knowledge of the Casimir-Polder 
potential itself. 

The experiment could be used, in addition, to test the 
validity of our model. To this aim, measurements should 
be performed in the near regime (say within the first 
ten wells) with a single isotope. In this case, the mea- 
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FIG. 8: (Color online) In yellow are displayed the regions 
of the (oy, Ay) plane excluded by experiments. The figure 
is taken from !4q|. The three superposed curves represent 
the experimental constraints theoretically calculated for the 
experiment FORCA-G. They correspond to the near regime, 
using a superposition between wells n = 4 and n = 6 (blue 
solid line, first from the left), the far regime for n = 40 (red 
solid line) and for n = 70 (black dashed line). 

sured energy differences would check the consistency of 
our atomic description as well as provide an estimation 
of the effective radius. The use of a single isotope makes 
the correction coming from the Yukawa potential negli- 
gible with respect to the Casimir-Polder term (see tables 
inilandHVl). 

Finally, the experimental setup can be used for a mea- 
surement of the Casimir-Polder potential around 5 fim: 
in this region, as shown in section |Vl the energy correc- 
tion due to the atom-field interaction is almost insensitive 
to the model chosen and the Yukawa interaction is much 
smaller than the quantum electrodynamical one. This 
measure could provide a new experimental observation of 
the Casimir-Polder potential with a relative uncertainty 
of less than one part in 10^. 

Nevertheless, we stress here that a precise knowledge of 
the Casimir-Polder standard potential requires an accu- 
rate description of the atomic and surface optical data. 
The details of the latter are unavailable at present, so 
the calculation in this paper have performed for perfectly 
conducting mirror. To complete our analysis it will be 
enlightening to conipare our results to the exact sphere- 
plate calculations [42l, '43'] . In this remarked in 
[43| , an appropriate description of the dielectric proper- 
ties of the sphere is needed to mimic the atomic optical 
response. 

VII. CONCLUSIONS 

In this paper we have discussed the modifications of 
the Wannier-Stark states in presence of a surface. As 



a first step we have considered the presence of the sur- 
face as a boundary condition of the time-independent 
Schrodinger equation obtaining in this way a new class 
of states. These states, even if asymptotically coincident 
with the ordinary Wannier-Stark states at large distances 
from the surface, significantly differ from them, both in 
energy and shape of the wavefunction, at the first few 
wells. 

We have then also taken into account the Casimir- 
Polder interaction between the atom and surface as a 
source of correction to the energy levels of the system. 
We have shown that these corrections diverge due to the 
behavior of the electrodynamical potential energy. 
In order to regularize this result, we have introduced an 
effective description of the atom as a probability den- 
sity distributed over a spherical volume. Our description 
leaves as free parameters both the radius of the sphere 
and the probability distribution. We have characterized 
the dependence of our results on both quantities. The 
validity of this model as well as the values of these pa- 
rameters remain to be investigated by experiments. 

In the second part of the paper we have studied the 
possibility of measuring a hypothetical Yukawa-type con- 
tribution to the gravitational potential at short distances. 
We have calculated the constraints that the experiment 
FORCA-G will be able to set on the (ay, Ay) plane. We 
have shown that the constraints set by the experiment are 
dominated by the experimental uncertainties and unaf- 
fected (to within those uncertainties) by the choice of the 
model for the regularization of the CP interaction. 

This work paves the way to the precise calculation of 
the energy levels in the experimental configuration of 
FORCA-G and other experiments that use atoms in op- 
tical dipole traps close to a surface P, 0| . To this aim, a 
precise knowledge of the optical data of the mirror and 
the atom is needed. This information will allow us to 
give a more detailed estimate of the accuracy of our re- 
sults, also based on the comparison with independent 
approaches to the regularization problem. Finally, the 
knowledge of the atomic wavefunctions constitutes the 
first ingredient for the description of the dynamics of the 
system, which is the subject of ongoing work. 
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